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ABSTRACT 

Context. Star formation depends strongly both on the local environment of galaxies, and on the internal dynamics of the interstellar 
medium. To disentangle the two effects, we obtained, in the framework of the AMIGA project, H„ and Gunn r photometric data for 
more than 200 spiral galaxies lying in very low-density regions of the local Universe. 

Aims. We characterise the H„ emission, tracing current star formation, of the 45 largest and less inclined galaxies observed for which 
we estimate the torques between the gas and the bulk of the optical matter We could subsequently study the H„ morphological aspect 
of these isolated spiral galaxies. 

Methods. Using Fourier analysis, we focus on the modes of the spiral arms and also on the strength of the bars, computing the torques 
between the gas and newly formed stars (H„) and the bulk of the optical matter (Gunn r). 

Results. We interpret the various bar/spiral morphologies observed in terms of the secular evolution experienced by galaxies in 
isolation. We also classify the different spatial distributions of star forming regions in barred galaxies. The observed frequency of 
particular patterns brings constraints on the lifetime of the various evolution phases. We propose an evolutive sequence accounting 
for the transitions between the different phases we could observe. 

Conclusions. Isolated galaxies appear not to be preferentially barred or unbarred. Through numerical simulations, trying to fit the 
Hj, distributions yields constraints on the star formation law, which is likely to differ from a genuine Schmidt law. In particular, it is 
probable that the relative velocity of the gas in the bar also needs to be taken into account. 

Key words, galaxies: general - galaxies: evolution - galaxies: kinematics and dynamics - galaxies: spiral - galaxies: statistics 

1. Introduction stellar bar (e.g.. Ide Vaucouleurs & de Vaucouleursl Il963l) . The 

gravitational torque of the bar makes the gas lose angular mo- 
The presence of a bar in a spiral galaxy is a striking feature and, mentum, driving it towards the centre and creating a central mass 
as such, IS one of the fundamental elements of the first morpho- concentration. The latter is able to perturb the elongated orbits 
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H , logical classification done by Hubble. The bar (including ovals) supporting the stellar bar, deflecti ng the stars p assing close to 
. . ^ , frequency determined in galaxies depends on the wavelength of ^^e centre , so weakening the bar (iPfenniger & N orman. 199C|. 
the images but is always higher than 65%. Some estimations However. iB^aud et al.1 (EoOsh have shown, using fufly self- 
based on NIR images, not aff^ected by extinction and tracing consistent simulations, that with gas parameters typical for nor- 
mainly the old population, reveal that_as many as 90% of the j^^j ^pj^^j^^ ^^e mass concentration i s not sufficient to fully di s- 
galaxies have bars dSeigar & JamesHl998|) . sol^g j^e bar, as was also claimed bv lShen & Sellv^ Jlool . 

Numerical simulations have establish ed that bars in gas-rich 
spiral galaxies are short-lived structures (iBournaud & CombesL 

l2O02h . At least two mechanisms have been proposed to weaken The second proposed mechanism to weaken a bar are the 

flie bars. The first one is the building of a large central mass con- bar torques themselves. The bar is a negative angular momen- 

centration, due to the gas inflow to the cen tre through the neg- tum wave, while the gas exerts a positive torque on the bar, due 

ative torques exerted on the gas by the bar (lAthanassoula et al.L to the balance between action and reaction forces. The angular 

I200I IBournaud et all l2005h . These torques are proportional to momentum lost by ttie gas is gained by the bar, which dissolves 

tiie phase shift between the gas and the stellar bar. It is well progressivel y For typical Sb-Sc galaxies the bar is destroyed in 

known fliat the gas is concentrated on the leading side of the about 2 Gyr ( IBournaud & Combesj |2002l) The observatio n of a 

high bar frequency from z ~ to z ~ 0.7 jEskridge et al.L l20()2t 

* Partially based on observations collected at the Centro Astronomico Up^^^ aH |200l I Sheth et al. |, HqoH) cannot flius be interpreted 

Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- ^ support the existence of robust, long- lived bars. Instead, this 

Planck Institut fur Astronomic and the Instituto de Astroffsica de supports the frequent renewal of bars. BerentzeiLet al. (2004|) 

Andalucfa (CSIC), as well as at the Observatorio de Sierra Nevada have shown that interactions can only form bars in gas-poor 

(OSN), operated by the lAA/CSIC. galaxies, which is not the case for most spiral galaxies. 
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Bar renewal can occur when the disks of spiral galaxies 
are replenished in cold gas through external accreti on, able 
to inc rease significantly the disk to bulge ratio. Bl ock et al.l 
(I2OO2I) suggest that external accretion of gas in the disk of 
spirals plays a fundamental role in explaining the high frac- 
tion of barred spirals and the observe d torque dis t ributio n. For 
instance, a recent study of M 33 by iBlock et al] ( |2004|) illus- 
trates that gas is feeding the galaxy through external accre- 
tion. Isolated, non-accreting spiral galaxies are very unlikely. To 
better understand the frequency of bars and their origin, mea- 
surements of the gravitat ional torques and bar forces in field 
galaxies have been made felock et aTl 120021: iLaurikainen et al.L 
[2004a. b; iButa et al.L |2004|) . but not for a sample of isolated 
galaxies. Such a sample is provided by the AMIGA project 
(Analysis of the interstellar Medium of Isolated GAlaxies), 
which constitutes a statistically significant template in the study 
of star formation and galaxy evolution in denser environments. 
We used as a start ing point the Catalogu e of Isolated Galaxies 
(CIG) compiled bv lKarachentseval ( 1 1 973h (see Sect.|2]i for which 
a careful re evaluation of degree of isolation is presented in 
IVerlev et aLl fcOOTa b). Previous AMIGA papers evaluate, refine 
and im prove the sample in diff'erent ways i ncluding revised po- 
sitions JEeo^^^^rfe^Montenei^ 2003|), optical characteri- 
sation ( Verdes-Montenegro et al.L 120051) . morphological revision 
dSulentic et al.L |2006|) . analysis of mid- and fa r-infrared prop- 
erties using data from the IRAS survey ( iLisenfe ld et al., 2007) 
and a study of the neutral (CO and H i) gas ( Espada et all 120051 : 
lEspadal [20061) . We also plan to study radio-continuum emission 
to determine current star formation rate and possible nuclear ac- 
tivity. The data are being released and periodically updated at 
|http : / /www . amiga . iaa.es 

In this paper we present a study of a carefully selected sub- 
sample of AMIGA, which is described in detail in Sect. [2] The 
observations and main reduction steps are summarised in Sect. [3] 
Then, we analyse, through Fourier analysis of gravitational po- 
tential and density, the intensity of the perturbations in each of 
the 45 galaxies (Sect.[4]i. Numerical simulations (Sect. [5) of stars 
and gas in isolated galaxies give us the theoretically expected 
star formation. Our results are discussed in Sect. [6] and the con- 
clusions are summarised in Sect. [7] 



2. Sample of isolated spiral galaxies 

As explained above, the AMIGA cata logue is based on th e 
Catalogue of Isolated Galaxies (CIG, iKarachentseval Il973l) . 
This catalogue is a compilation of 1051 objects with appar- 
ent B magnitu d e brig hter than 15.7 and declination > -3°. 
iKarachentseval ( Il973h visually inspected the Palomar Sky 
Survey prints, trying to identify those galaxies in the Catalogue 
of Ga laxies and Clusters of Galaxies (CGCG, IZwickv et al.L 
[T968h which have no near neighbours. An extensive study of the 
isolation of these galaxies is presented in V erlev et al . (2007a b), 
as a part of the AMIGA project, and discussed in Sect. [6] 

One of the aspects of our project was to study the Hq, emis- 
sion of isolated spirals in its different aspects. With this aim, and 
in order to avoid well known biases due to flux- or magnitude- 
limited samples and to end with a complete and homogeneous 
sample, we kept all the galaxies in a volume-limited sample, i.e. 
with observed recession heliocentric velocities Vr: 

1500 km s ' < Vr < 5000 km s ' . 

It represents about one fourth (251 galaxies) of the whole 
CIG. Among them, 27 were early-type galaxies according to our 



new morphological classification dSulentic et al.L l2006h . hence 
rejected from the sample. We finally ended up with 224 galax- 
ies with morphological types SO/a or later. While being a rather 
big sample allowing statistical studies of the properties of iso- 
lated spiral galaxies in the local Universe, its size still allows to 
perform H,, imaging for all the galaxies. These 224 galaxies con- 
stitute the main subsample for AMIGA star formation studies. 

For the here presented study two further requirements had to 
be imposed: 

1. Major axis a> 1' to have a sufficient spatial resolution; 

2. Inclination ; < 50° in order to obtain a sufficiently accurate 
deprojection. 

From our subsample, 48 galaxies follow these two criteria, 
of which 45 were observed and which constitute the subject of 
study in this paper. Table ?? reports information about the main 
characteristics of the selected galaxies, along with some techni- 
cal details. The position angle (PA), inclination / and rotation di- 
rection (assuming trailing spiral arms, the rotation will be set as 
positive if counter-clockwise) listed are those used to calculate 
the torques (24 counter-clockwise galaxies and 19 clockwise). 

3. Observations and data reduction 

The 224 galaxies comprising our sample for star formation stud- 
ies, including the 45 studied in this paper, have been imaged 
in H„ and Gunn r filters, in order to trace respectively the re- 
gions ionised by newly born stars (Hii regions) and the older 
stellar component. The observations and data reduction will be 
described in more detail in a separate paper (Bergond et al., in 
prep.), and are just summarised here. For each galaxy, we se- 
lected the more appropriate red-shifted H„ filter (narrow-band 
filter, typically 50 A), based on its observed recession velocity 
and we used the Gunn r filter (broad-band filter, typically 870 A) 
to subtract the continuum contribution affecting the narrow -band 
filter The Gunn r broad-band images trace the stellar component 
of the galaxy, and the H„ images show the young stars born from 
the gas (Hii regions). The individual exposure times were 300 
seconds for the Gunn r filter and 1200 seconds in H„. We ap- 
plied small shifts between two successive exposures to be able 
to treat bad pixels on the chip. We applied the usual corrections 
(bias subtraction and division by twiUght flat fields) to remove 
the instrumental signature. 

The removal of the continuum contribution to the flux in the 
images taken through the narrow- band Hq. filter is the most deli- 
cate task of the reduction process. iJames et"al1(l2004h tested var- 
ious methods: numerical integration to find the ratio of the fil- 
ter profile integrals, photometry of standard spectrophotometric 
stars through the pair of filters to find the scaling factor, use of 
foreground stars in the narrow-band and continuum images of 
each galaxy. The two last methods gave the most consistent and 
accurate results. The latter method has two advantages: it takes 
into account any changes in the sky transparency between the 
two images and the possibility to use several stars improves the 
statistics and leads to the most accurate removal. Hence, we ap- 
plied this last method to find the scaling factor to be applied 
to the images used for continuum subtraction. As the H„ wave- 
length is comprised in the Gunn r filter, a significant H,, emis- 
sion would perturb the continuum removal. The uncertainties in- 
volved can be estimated in various ways. In a first approxima- 
tion, with the hypothesis of negligible line emission (assuming 
a flat spectrum), the ratio of the width of the filters shows that 
the uncertainties amounted to about 6% (using typical width of 
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50 and 870 A). Photometry of the total emission on the Gunn r 
and Ha-continuum images confirms that the H„ emission gen- 
erally represents a small fraction (typically 3%) of the Gunn r 
e mission. This latter v alu e agrees with the esti mations reported 
bv lJames et al. I (12004 and lKnapen etall (|2004 . 



4. Image analysis 

In order to calculate the gravitational potential, density and 
torques in the plane of the disks of the galaxies, we processed 
the H„ and Gunn r corrected images as follows: the r band im- 
age (not scaled to H„) and the H„ continuum subtracted image 
were used. The Gunn r images were used in order to estimate 
the gravitational potential of the galaxies. This step is gener- 
ally done with near-infrared images which trace more closely 
the gravitational potential (absence of extinction); unfortunately 
we do not possess these images for all the galaxies in our sample. 
Nevertheless, CIG 100 4 (NGC7479) is in the OSUBGS sample 
dEskridge et al.L 2002h an d its near-infrared potential had been 
estimated bv iBlock et al.l (|2P02). The potentials obtained from 
the Gunn r and near-infrared images are very similar (only the 
potential due to the optical spiral arm shows a slight enhance- 
ment in the Gunn r with respect to the near-infrared estimation) 
and hence the Gunn r image is a rather good approximation in 
order to derive the gravitational potential of a galaxy, in the ab- 
sence of near-infrared information. We defined the centre as the 
maximum in luminosity near the geometrical centre. The cen- 
tre is the same in the r and H„ image to avoid artificial torques 
between the two components. Then, we cut 512 x 512 pixels^ 
subsections, as Fast Fourier Transform (FFT) will be applied on 
these images. Foreground galactic stars were removed in order to 
avoid a contamination in the density and potential derived from 
the luminosity of the pixels. The galaxy was deprojected to ob- 
tain the gravitational potential in the disk plane. The deprojec- 
tion assumes that the outer disk of galaxies are circular and that 
the bulge s are as flat as the disks. For galaxies inclined by less 
than 70°. lBlock et alj ( l2002h found that the observational uncer- 
tainties on the incUnation angles of the galaxies affect barely the 
estimation of the torques. Our conservative limit includes galax- 
ies inclined less than 50° and assures us that the uncertainties in 
the deprojection do not affect our results. 

Then, to determine the bar/arm force at each radius, we made 
use o f the method developed by [Garcia-Burillo et al. (1993^ 
l2005b to compute the gravitational potential, applying it to our 
red images, supposing a constant stellar mass-to-light {MIL) 
ratio. A Fourier-component analysis of the potential was per- 
formed, using the m - 2 component to get the tangential forces 



dBlock et 



smgtnem 
^ l2002h 



and the m - component to get the ax- 
isymmetric forces. Their quotient gives us a measurement of the 
bar/arm force. 

We estimated the average torque depending on the radius: 
using the gravitational forces and the young stars born from the 
gas (Hq,), we obtained the phase shift between gaseous arms 
and the potential well which creates the torques. Generally, 
due to the bar gravity torques, the gas flows into the centre 
(negative torque) from the corotation radius to the inner Linblad 
resonance (ILR) and flows outwards (positive torque) when 
outside the corotation radius, until the outer Lindblad resonance 
(OLR). This can be elTectively observed for the galaxies in our 
sample. Finally, we estimated the angular momentum transfer 
and therefore the evolution time for bars in isolated galaxies 
(see subsection l6.4.2l i. 



The gravitational potential was decomposed as: 

0(r, 6) = (I'o(r) + ^ (D„,(r) cos(»z6' - (^,„) 

111 

where the strength of the m-Fourier component is 
Qiiiif) = m<[>„,/r|Fo( r)|, and its global streng th over the 
disk: max. (Qmir)) (e.g.. (Combes & Sandersl[T98l . 

The (disk) surface density was decomposed as: 

ju(r, (p) = ^iQ{r) + ^ a,„(r) cos{}n(f> - <f>,„(r)) 

m 

where the normalised strength of the Fourier component m is 
A„,{r) = a,„/fio(r). The density Fourier analysis is complemen- 
tary to that of the potential, which is less noisy but more global, 
and depending on the assumed 3D distribution of the mass. 
The maximal torque at a given radius was defined by: 



Qt{R) 



77 max 



(R) 



Fq(R) 



1 / ao(«,9) \ 
R\ ae In 

dOo(«) 
dR 



where Fj™(R) represents the maximum amplitude of the 
tangential force at radius R, and Fo{R) is the mean axisymmetric 
radial force inferred from the m - component of the gravita- 
tional potential. 

As an example, Figs.[T]&|2]display the images and graphics 
of the potential, density and torques corresponding to CIG 30 & 
CIG 1004, respectively. The first row shows the Gunn r broad- 
band filter (left) and H„ narrow-band filter (right) deprojected 
images, from which the potential, density and torques have been 
estimated. The second row corresponds to the amplitude (left) 
and the phase (right) of the Fourier components of the poten- 
tial. The third row corresponds to the amplitude (left) and phase 
(right) of the Fourier components of the density. For the graph- 
ics of the potential and the density, the legend is the following: 
dashed line represents m = 1; full line represents the m - 2 
component (and the sum also for the amplitude of the poten- 
tial as it is always above and cannot be confused with the oth- 
ers); dot-dash-dot-dash represents m - 3; dotted line represents 
m - 4. The torques as a function of the radius of the galaxy 
is presented in the bottom left panel. The corotation occurs at 
the radius defined by the change of sign for dL/L. And dL/L 
shows that in one rotation, the gas loses a non negligible frac- 
tion of its angular momentum, which can bring constraints on 
the life time of the morphological patterns of the galaxy. The 
legend in the bottom right corner gathers the CIG number, the 
observed recession velocity (in km s"'), the optic al diameter D25 
in 0, the blue luminosity (logarithm, in L©; see iLisenfeld et all 
l2007h . their Table 3), the morphological type (from NED) and 
the group to which the galaxy belongs (see Sect.|6]l. The Gunn r 
and H„ images, along with the graphics of the potential, density 
and torques of the full sample of the 45 galaxies are also avail- 
able at http : //www . aitiiga. iaa. es/publications[ 



4 



S. Verley et al.: Star formation in isolated AMIGA galaxies: dynamical influence of bars 



Table 1. Initial conditions for the numerical simulations. A^p is 
the number of particles considered, rb/d/h and Mb/d/h the galaxy 
bulge/disk/halo radii and masses, and /ei the collision elasticity 
factor. 



Run 






Mb 


'■d 


Mi 






/el 


id. 


xlO^ 


kpc 


(1) 


kpc 


(1) 


kpc 


(1) 




RunO 


100 


1.1 


2.5 


4.4 


8.0 


16.0 


7.2 


0.65 


Run 1 


100 


1.1 


6.8 


4.4 


8.0 


16.0 


11.8 


0.65 


Run 2 


100 


1.1 


6.8 


4.4 


8.0 


16.0 


11.8 


0.85 


Run 3 


400 


1.1 


6.8 


4.4 


8.0 


16.0 


11.8 


0.65 


Run 4 


400 


1.1 


6.8 


4.4 


8.0 


16.0 


11.8 


0.85 



'"Masses are in lO'^M©- 



5. Numerical simulations 

5.1. Methods 

In order to understand the observed H„ distributions, and the dif- 
ferent phases identified, we performed A-body simulations with 
stars and gas, including star formation. Since we want to explore 
many physical parameters, we chose to carry out 2D simulations, 
which should capture the essential of the bar evolution, and lo- 
cation of star formation in these isolated galaxies. The 3D com- 
ponents, bulge and dark matter halo, are therefore considered 
as rigid and spherical potentials, in which the disk component 
evolves. The rigid or live character of the spheroidal components 
would make a significant effect, by exchanging angular momen- 
tum, only if their mass is relatively significant. But in the inner 
parts of galaxies, the dark matter halo is not dominant, and thus 
this phenomenon would only induce small changes. The same is 
true for the spheroidal bulges, which are relatively small in mas, 
in these mostly late-type galaxies. Self-gravity is only included 
for the disk (gas -i- stars). 2D A-body simulations were carried 
out using the FFT algorithm to solve the Poisson equation, with 
a Cartesian grid, varying from 256 x 256 to 512 x 512 (useful 
grid, free of periodic images). Two spatial resolutions were se- 
lected to appreciate its influence on the star formation physics. 
The cell size is then from 62.5 to 125 pc, and the total size of 
the grid is 32 kpc. The softening length of the gravity has the 
characteristic scale of the cell (62.5 to 1 25 pc). More details o n 
the numerical techniques can be found in lCombes et alj(ll990h . 

The stellar component is represented by 100000 or 400000 
particles, and the gas component by 40000 and 160000 for the 
low and high resolutions adopted, respectively. 

The bulge is modelled as a rigid spherical potential with 
Plummer shape: 

Ob(r) = —= 

Mb and are the mass and characteristic radius of the bulge. 

The stellar disk is initially a Kuzmin-Toomre disk of surface 
density: 

truncated at 15 kpc, with a mass M^. It is initially quite cold, with 
a Toomre Qt parameter of 1 . The halo is also a Plummer sphere, 
with mass Mh and characteristic radius r^. The time steps are 
0.5 and 1 Myr. The initial conditions of the runs described here 
are given in Table \T\ The particle plots of the stars and gas as 
well as the recently formed stars are depicted in Figs. |3]&|4] and 
Figs.|5]&|6] for Run 1 and Run 4, respectively. 



The gas is treated as a self-gravitating component in the A- 
body simulation, and its dissipation is treated by a sticky particle 
code, as in Combes & Gerin ( 1985). The initial gas-to-total mass 
ratio (Fgas) in the disk ranges between 6 and 14%, since the star 
formation in the simulation is capable of reducing Fgas to a fi- 
nal lower value. The mass of one gas particle therefore varies 
between 8 x 10"* and 3 x IO^Mq. 

The initial distribution of gas in the model is an exponen- 
tial disk, truncated at 15 kpc, and with a characteristic radial 
scale of 6 kpc. Initially, its velocity dispersion corresponds to a 
Toomre Qj parameter of 1. The gas clouds are subject to inelas- 
tic collisions, with a collision cell size between 60 and 120 pc 
(region where particles are selected to possibly collide). This 
corresponds to a lower limit for the average mean free path of 
clouds between two collisions. The collisions are considered ev- 
ery 5 to 10 Myr In a collision, the sign of the relative cloud 
velocities is reversed and the absolute values are reduced: rela- 
tive velocities after the collision are only /d times their original 
value, the elasticity factor /ei being between 0.65 and 0.85, as 
indicated in Table [T] The dissipation rate is controlled by this 
factor. All gas particles have the same mass. 

Star formation is taken into account following a generalised 
Schmidt law: the star formation rate is proportional to the vol- 
ume density to the power « = 1 .2, provided that the gas volumic 
density is larger than pgas = 1 H-atom cm"-', i.e. the rate of gas 
mass transformed into stars is dm oc pl'.^^dt. To compute this rate, 
at regular intervals of df = 5-10 Myr, the gas density is aver- 
aged in each cell, and the probability of the gas particles being 
transformed into stars is computed by: 

P = d;7z/Meell 

for all particles in this cell, of mass Mceii. Each new star formed 
has exactly the same mass as each gas particle, about 3 times 
smaller than any old stellar particle. This simple scheme corre- 
sponds to an instantaneous recycling of matter, since the contin- 
uous mass-loss from recently formed stars is not followed. The 
rate of star formation is normalised so that in unperturbed runs 
(without galaxy interaction, galaxies are quiescently and regu- 
larly forming stars), the timescale for consumption of half of the 
gas mass is of the order of 2 Gyr (SFR ~ 1-2 MQ/yr). At each 
star formation event, the neighbouring gas particles are given a 
small extra velocity dispersion of the order of ~ 10 km s ' . 

To reproduce a model equivalent of the H„ maps, we plot 
only the recently born stars, still able to ionise their surround- 
ings. In order to have more statistics, we choose to plot all new 
stars of age less than 200 Myr (Figs.|4]&|6ll. 

5.2. Results 

We tried several runs corresponding to the various morphologi- 
cal types of the Hubble sequence, varying the bulge to disk mass 
ratio, the halo to disk mass ratio and the gas content. This had the 
consequence to vary the dynamical stability of the disk: the bar 
appears within different time-scales, and with mainly two differ- 
ent morphologies for, respectively, early-types and late-types. In 
the case of late-types, the bar develops without any resonant ring, 
and quickly the gas infalls to the center, and destroys the bar 
The H„ morphology has then a bar shape during a short period, 
and quickly peaks at the center. For early-types, the high central 
mass concentration allows the formation of inner Linblad res- 
onances, and gaseous resonant rings are developped inside the 
bar Also the life-time of the bar is long enough to develop a 
resonant ring at the 4: 1 resonance near corotation, surrounding 
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Fig.l.IC35 (CIG0030). 

V =4586kms-' 
Optical diameter: 1.0 arcmin. 
Blue luminosity: 9.82 
Morphology: Scd: 
Group: E 

For the potential, Q\, = Q„„ the normalised amplitudes of the 
m Fourier modes. For the graphics of the potential and the 
density, the dashed line represents m = 1 ; full line represents 
the m = 2 component (and the sum also for the amplitude 
of the potential as it is always above and cannot be confused 
with the others); dot-dash-dot-dash represents m = 3; dotted 
line represents m = 4. 
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Fig. 2. NGC 7479 (CIG 1004). 

V = 2376kms-' 
Optical diameter: 4. 1 arcmin. 
Blue luminosity: 10.40 
Morphology: SB(s)c; LINER Sy2 
Group: G 

For the potential, Qb = Q„„ the normalised amplitudes of the 
m Fourier modes. For the graphics of the potential and the 
density, the dashed line represents m = 1 ; full line represents 
the m = 2 component (and the sum also for the amplitude 
of the potential as it is always above and cannot be confused 
with the others); dot-dash-dot-dash represents m = 3; dotted 
line represents m = 4. 



S. Verley et al.: Star formation in isolated AMIGA galaxies: dynamical influence of bars 



7 




Fig. 3. Particle plots of the stars (top) and gas (bottom) at the 
epoch r = 1.1 Gyr (left) and T = 1.5 Gyr (right), in the simula- 
tion of Run 1. The scale of the square is 64 kpc. 



Fig. 5. Particle plots of the stars (top) and gas (bottom) at the 
epoch T - 0.6 Gyr (left) and T - \ .Q Gyr (right), in the simula- 
tion of Run 4. The scale of the square is 64 kpc. 




Fig. 4. Plots of recently formed stars (less than 200 Myr ago), 
at the epoch T - l.l Gyr (left) and T - 1.5 Gyr (right), in the 
simulation of Run 1. The scale of the square is 32 kpc. These 
plots should give an idea of the Hq, map. 




Fig. 6. Plots of recently formed stars (less than 200 Myr ago), 
at the epoch T = 0.6 Gyr (left) and T = 1.0 Gyr (right), in the 
simulation of Run 4. The scale of the square is 32 kpc. Note the 
orientation of the young star component, which keeps the phase 
shift of the gas component during some time, and is leading with 
respect to the old stellar bar of Fig.|5]by about 3 degrees. 



the bar (Fig. |3]l. However the gas density there is not enough to 
produce substantial star formation (Fig.|4]i. 

Whatever the morphological type of the galaxy, the mod- 
els always show a first stage where the main star formation is 
aligned with the bar, and a second stage, where the star forming 
regions concentrate towards the nucleus. Variations of the degree 
of dissipation of the gas component, through the variation of col- 
lision elasticity /ei, does not change the picture, as shown by the 
comparison between Run 1 and Run 4 in Figs. |4]and|6] The cor- 
responding H„ maps that can then be predicted are either of the 
group G, as defined in the next section, or the group H. But it is 
not possible to obtain the most frequently observed morphology 
of the group E. In other words, the equivalent Hq, maps are pre- 
dicted to correlate tightly to the gaseous maps. This is expected 
from the type of recipe adopted for the star formation rate, i.e. a 
Schmidt law, with a density threshold. We defer the discussion 
of these results to Sect. 16.4.21 



6. Discussion 

Our sample being 94% complete (we are only missing data for 3 
galaxies over the 48 we had selected, which could not have been 
observed), we can describe statistically significant characteris- 
tics for a sample of isolated galaxies. 

6.1. Isolation of the galaxies 

IVerlev et~al] (l2007allbh lead an extensive study to revise the de- 
gree of isolation of 950 CIG galaxies. Within a projected ra- 
dius of 0.5 Mpc, all the neighbours that could exert an influ- 
ence on the evolution of the CIG galaxies were catalogued. 
Parameters quantifying the degree of isolation of each CIG 
galaxy were defined. A particularly relevant parameter, Q, mea- 
sures the strength of the tidal forces affecting the central galaxy. 
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normalised to the inner binding forces of the galaxy itself. In the 
subsample of 45 galaxies defined in the present study, three of 
the galaxies present tidal forces amounting to more than 10% 
{Q > -1) of their inner binding forces. This means that their 
evolution could have been seriously affected by nearby compan- 
ions. These three galaxies are CIGs 80, 85 and 1001. 

Also, other eleven galaxies present values of the external 
forces amounting to more than 1% of their internal binding 
forces {Q > -2), which is the limit from which the environ- 
ment can partially in fluence the evolution of the central galaxy 
anassou lal ll984 . These galaxies are CIGs 59, 84, 96, 176, 
376, 382, 512, 645, 808, 812 and 1004. Nevertheless, these 
galaxies are isolated from similar size neighbours (wi t hin a fac- 
tor 4 in size, as originally defined bv iKarachentseval Il973[) . A 
statistical rejection of the projected distant neighbours (with a 
difference in recession velocity of almost 20 000 km s"') reveals 
that all but one (CIG 176) can, indeed, be galaxies whose evolu- 
tion is dominated by their internal properties. More data on the 
recession velocities of the neighbours are needed in order to ac- 
cess a definitive picture of the environment of each CIG galaxy. 

6.2. Percentage of bars 

To know the percentage of barred galaxies in a sample of isolated 
galaxies, we visually classified the 45 galaxies in three subsam- 
ples, as a function of the presence or absence of bar: 

Bar 27 galaxies - CIGs 30, 53, 66, 96, 116,1 76, 188, 250, 267, 
376, 382, 512, 645, 652, 660, 661, 712, 754, 808, 840, 854, 
862, 931, 935, 1001, 1004, 1039. 

No bar 15 galaxies - CIGs 50, 59, 80, 84, 85, 217, 281, 291, 
359, 700, 750, 812, 875, 924, 992. 

Intermediate 3 galaxies - CIGs 68, 575, 744. 

The barred galaxies represent 60% of our sample (the bar can 
be seen either in Gunn r or in H„). As we are using optical im- 
ages, this result is consistent with the studies presented in Sect.[T] 
The unbarred galaxies represent 33% of the sample. Three galax- 
ies have an intermediate stage bar, similar to the "SAB" kind, 
they represent 7% of the total. As the morphological recognition 
is highly dependent both on the material inspected and on the 
human factor, some of the galaxies that we classified as "barred" 
would have shifted into t he "intermediate" cate gory following 
the classification done by Ide VaucouleursI (Il963h . Nevertheless, 
we can infer that the isolated galaxies span the whole range of 
bar morphologies, in quantities similar to the galaxies in denser 
environments. Isolated galaxies are not preferentially barred or 
unbarred gala xies. T his result is marginally in contradiction with 
the study by Yarela et a l - (2004), who estimate that bars are twice 
more frequent in perturbed galaxies compared to isolated g alax- 
ies, especially for early-types. Also lEknegreen et all (1 19901) find 
more bars in a sample of binary galaxies, and also more early- 
types. 

6.3. Phase shift between gas and stellar components 

We now focus on the existence of shift angle between the lo- 
cation of the gas and the location of the older stellar compo- 
nent, which has already been noted by iMartin & Friedlil (Il997h 
in a sample of eleven galaxies (for other exarnples, s ee also 
Ide Vaucouleurs & de VaucouleursI 119631 : iPhilUpslI 19961) . It was 
sometimes difficult to define where exactly the arms begin be- 
cause the H„ could be very clumpy, but frequently, we saw bright 
H„ knots at the end of the bars, even when there was no emission 



in the bar Consequently, we could use the star formation spots 
at the end of the bars, to define the starting point of the spiral 
arms. The H„ emission was always leading with respect to the 
bar in the Gunn r image. The most evident cases of bar shifts 
between the Gunn r and the H„ images are listed hereafter (in 
parentheses: an estimation of the shift angle in degrees): CIG 30 
(30°); CIG 53 (10°); CIG 66 (5°); CIG 96 (10°); CIG 176 (10°); 
CIG 376 (30°); CIG 512 (10°); CIG 840 (15°); CIG 1004 (5°). 
Three cases are illustrated in Fig.|7] where the orientations of the 
H„ and stellar bars have been marked by dashed black lines. 

6.4. Evolutive sequence 
6.4.1. Classification 

As some characteristics were frequently found among the galax- 
ies of our sample, we chose to group the galaxies presenting sim- 
ilar features. In decreasing frequency order, we defined the three 
main groups: 

Group E 19 galaxies - CIGs 30, 53, 66, 80, 96, 176, 217, 281, 
291, 376, 382, 660, 812, 840, 862, 931, 992, 1001, 1039. 
The principal features of this group are the following: a 
strong central peak in the Ha emission; no Hq, emission in 
the bar (for the barred galaxies in the Gunn r image); bright 
H„ knots at the end of the bar/beginning of the spiral arms; 
Ha emission along the spiral arms, generally clumpy. 

Group F 9 galaxies - CIGs 50, 59, 84, 188, 267, 661, 808, 924, 
935. 

This group is constituted by galaxies with less gas, having 
a smoother morphology. These galaxies do not present any 
central emission spot in H„. 
Group G 8 galaxies - CIGs 116, 250, 645, 700, 712, 754, 854, 
1004. 

This group gathers galaxies presenting Ha emission in the 
bar CIG 700 presents a very faint emission, but a bar can be 
distinguished. 



Five of the galaxies of our sample did not fit in any of the main 
groups defined above, while four others presented characteristics 
mixing features from two of the groups: 

Group H 3 galaxies - CIGs 68, 359, 875. 

This group is mostly constituted by early types galaxies, with 
very few emission in H„ and when it occurs, only in the cen- 
tre. 

Group Irr 2 galaxies - CIGs 85, 575. 

This group gathers very irregular galaxies which do not fit 
particularly in any of the previous groups defined above. 

Group EG 3 galaxies - CIGs 512, 652, 744. 

This group is a transition between the E and G groups pre- 
sented above: the galaxies exhibit fragments of bars. The 
central emission of CIG 744 is very clumpy but seems to 
follow the shape of a bar, this is why we included it in this 
group. 

Group EE 1 galaxy - CIG 750. 

CIG 750 does not have a clear central peak in but 
presents all the other features of the galaxies in the group 
E: clumpy H„ emission along the spiral arms. We call this 
group EE as a result of the mixed features presented by 
CIG 750. 
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Fig. 7. Examples of shift angles between the H„ and stellar bars. The dash line represents the positon angle of the perturbation 
found in the image analysis. Left: Deprojected images of CIG 30 in H„-continuum (upper panel) and Gunn r (lower panel). Middle: 
Deprojected images of CIG 376 in H„-continuum (upper panel) and Gunn r (lower panel). Right: Deprojected images of CIG 840 
in Ha-continuum (upper panel) and Gunn r (lower panel). 



6.4.2. Interpretation 

Many A^-body simulations of galaxies with gas have re- 
vealed that bars are not long- liv ed in late-type galax- 



ies (e.g. iFriedli & Ben"3. 119931: iBerentzen et al.L Il998t 



iBournaud & Combesl |20M) . iBoumaud et alj ( l2005h have 
shown that the bar dissolution is due to the gas inflow towards 
the centre, and that bars can be renewed by external gas 
accretion. Our goal here is to get more insight in the various 
dynamical processes that could be involved in the interpretation 
of the different observed morphologies as different stages of an 
evolutive sequence. In our sample, we can identify various main 
steps occurring during the formation and evolution of the bars. 

By gravitational instability, a galaxy accretes gas from the 
inter-galactic medium which makes it unstable for bar formation. 
The bar creates a torque which drives the gas inflow towards the 
centre. This phase corresponds to our identified G phase (see, 
for instance, Fig.|2]for a typical galaxy belonging to this group). 

The second step is a transition between the G and E phases: 
the gas inflows towards the centre while a ring is slowly forming 
at the resonance (pseudo-ring due to the winding of the spiral 
arms, at the ILR). 

In a third step, the gas is progressively depopulated from the 
bar, and accumulates first in the very centre of the galaxy (or a 
very small nuclear ring, at ILR), and also at the ultra harmonic 
resonance (UHR), near the corotation. The gas there is quite sta- 



ble, with very small relative velocity with respect to the bar pat- 
tern, and the star formation could be quite efficient in these re- 
gions. This corresponds to our identified more frequent phase, 
the E phase, where galaxies could spend more than 40% of their 
life-span (see, for instance, Fig.[T]for a typical galaxy belonging 
to this group). 

Since the gas infall destroys the bar, the latter becomes pro- 
gressively weaker and weaker: the F phase is reached. The stars 
in the centre become an old population, contributing to increase 
the bulge mass. Without more gas fuelling, the H„ spot in the 
centre is fading away in 10^ years (OB stars). The frequency of 
the F phase means that a galaxy spends typically 20% of its life 
in this stage, where the bar is weakened or destroyed. 

The various classes proposed in the previous section were 
based on the H„ map morphologies. By plotting images of the 
recent stars formed in the simulations, we have tried to predict 
those H„ morphologies in Sect. |5] However, the two most fre- 
quent predicted morphologies are either corresponding to group 
G (Hq, emission aligned along the bar), or H, when the emission 
is concentrated in the nucleus. The configuration E, where the 
star formation avoids the bar, is not reproduced in the simula- 
tion, in contradiction to our observations where the E phase is 
about twice more frequent than the G phase. 

We think that this feature comes mainly from the star forma- 
tion recipe adopted in the simulation. The efficiency of the star 
formation as a function of the gas density is still a challenging 
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Fig. 8. Distribution of total non-axisymmetric strengths in the 
CIG sample. 
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Fig. 9. Distribution of total non-axisymmetric strengths in the 
OSU sample. 



issue. The predicted gas density appears to reproduce the ob- 
servations. Our simulations reveal the formation of rings at the 
UHR, but this is always accompanied with a large gas density 
along the bar 

In most of the barred galaxies, the gas is observed in the bar, 
through the dust lane or thr ough molecular line emission, like 
CO (e.g. iHelfer et all l2003h , but little or no star formation and 
Ho- is observed. With the Schmidt law used in the present simu- 
lations (SFR oc Pgj,';), it is expected to find recent star formation 
in the bar, while in the actual galaxies, t his is avoided, maybe 
through too large relative velocities (e.g. iRevnaud & Downei, 
Il998h . Therefore to remove the discrepancy between Hq, mor- 
phologies and predictions, it would be necessary to change the 
expression of the star formation rate: it should not only depend 
on pg but should take into account the relative velocity of the gas 
with respect to the bar. This will be presented in a future work. 



The absence of star formation and H„ emission in the bar, 
while the gas density is abundant, is a frequently observed phe- 
nomenon. We do not have galaxies in com mon with the BIMA 
Survey Of Neai'by Galaxies (BIMA SONG. lHeiferetan.l2002h . 
But already, the BIMA SONG brings us some clues about the 
gas density. The E phase, although also present there, is less fre- 
quent, relative to the G phase, as we could expect. There is more 
CO emission in the bars of their galaxy sample (NGC 2903, 
3627, 4535, 5457) compared to our Hq, sample. The H„ emis- 
sion is obviously a more efficient tracer of star formation than of 
the total amount of gas available. 

6.4.3. Comparison with otiier samples 

Other groups have also imaged and classified spiral galaxies 
as a f unction of their H„ morphology. Hameed & Devereu^ 
(Il999h . for instance, concentrated their work on 27 early-type 
spiral galaxies (Sa-Sab) in the nearby Universe. As our sam- 
ple involves mainly Sb-Sc galaxies, the two samples are some- 
how complementary, although it has to be noted that their 
galaxies are selected without any regard to isolation crite- 
rion. We can group their galaxies according to our classifi- 
cation. The E group would gather 11 galaxies: NGC 1350, 
NGC 1371, NGC 1398, NGC 1433, NGC 1515, NGC 1617, 
NGC2273, NGC3169, NGC5156, NGC7213 and 1108-48. 
Four galaxies would be classified in the G group: NGC 986, 
NGC 5188, NGC 7552 and NGC 7582. The central emission 
characteristic of the H group is seen in 5 galaxies: NGC 1022, 
NGC 1482, NGC3471, NGC3885, NGC5728. The galaxies 
NGC 972, NGC 5915 would be classified in the transitional EG 
group while it would not be possible to classify five galaxies 
(NGC 660, NGC 2146, NGC 3717, NGC 6810, NGC 7172) gen- 
erally due to their too high inclination angle which makes dif- 
ficult to disentangle the emission of the disk from the emission 
specific from the bar itself. 

The comparison between the two samples shows that three 
groups present the same percentages of galaxies: the percentage 
of group E galaxies is almost the same in the two samples (41 
and 42% in their and our sample, respectively). The same per- 
centage of galaxies is also found for the G group (15 and 18%, 
respectively) and for the EG group (about 7% in both samples). 
Additionally, all their ENER (Extended nuclear emission-line 
region) galaxies fall into our E group. Nevertheless, some 
differences are also observed. Our F group is not represented 
in their sample, which lacks totally the galaxies without central 
emission. The other notable difference is that their sample 
presents a proportion 3 times higher of H group galaxies. This 
can easily be explained because all their galaxies are Sa-Sab, 
hence more early-type (with less gas and very few emission) 
than the bulk of our sample. 

Other works tar geted the inner H„ emission of galaxies: 
iKnapen et al.l (l2006l) focussed their work on the central regions 
(within 2 pc) of 73 spiral galaxies, that can have relatively close 
companions, which is the case for 26 galaxies (36%) in their 
sample. We classified the 73 galaxies following our nomencla- 
ture: 45 galaxies (62%) are in the E group, 2 galaxies (3%) in 
the F group, 4 galaxies (5%) in the G group, 9 galaxies (12%) 
in the H group, 10 galaxies (14%) in the EG group, 1 galaxy 
(1%) in the EF group and 2 galaxies could not be classified (Irr 
group), partly due to their to high inclination. The restriction of 
their sample only to the 26 galaxies wich possess a physically 
linked companion does not change significantly the percentages 
of the galaxies in the groups (61% in E, 4% in F, 0% in G, 19% 
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in H, 12% in EG, 0% inE F and 1% in Irr). This is not surpris- 
ing as iKnapen et al.l ( l2006l) already noticed that the presence of 
a close companion does not produce any particular strong trend 
in the H„ morphology. 

The very low percentage of F galax ies with respe c t to o ur 
work is due to a selection criterion in the iKnapen et alj (l2006l) 's 
sample as they selected galaxies with some prior evidence for Hq, 
structure in their central regions. The same reason enhanced the 
proportion of group E galaxies because they selected galaxies 
with already known nuclear or inner rings. Their sample is not 
complete and the comparison and interpretation with our work 
should be done with care. With respect to the isolation, the cri- 
terion they used can not ensure that the galaxies without "very 
close companion" would be isolated following Karachentseva's 
criterion (even if it seems that the presence of a companion does 
not drastically repercuss on the Ha features). As a result we 
cannot compare directly our isolated galaxies with their sample. 
Nevertheless, the lHameed & Devereuxl(ll999l) and lKnapen et al.l 
( |2006|) studies show that the groups we defined are also seen in 
samples that might not be isolated, so that the evolutive sequence 
could also be applied to galaxies marginally in interaction. For 
more detailed comparisons, samples better defined in term of 
isolation (or interaction) and better statistics would be needed. 

6.4.4. Bar frequency 

We have computed the total strength of bar and spiral waves in 
our galaxy sample by the method outlined in Sect. |4] We have 
adopted an exponential scaleheight equal to 1/12 of the radial 
scale length, inde pendent of Hubble ty pe, consi stent with previ- 
ous studies (e.g.. lOuillen et all [T994I) although Ide Grij si (1 19981) 
found a type-dependent relationship and thicker disks. Varying 
this parameter will slightly change the absolute bar strengths in 
Figs. |8]&|9] but will not have much impact on the relative val- 
ues. The maximum of the parameter Qj, dominated by the bar 
strength and its harmonics, is used to build the histogram of the 
bar frequency for isolated galaxies (Fig.|8]l. It is then possible to 
compare with the same histogram (see Fig.|9]l build for a general 
sample, the OSU sample (Eskridge et al., 2002). This sample of 
galaxies, selected without regard to any isolation criterion, ap- 
pears to give quite similar results. Our histogram here is more 
uncertain and noisy, because of the much smaller sample size. 
However, it is possible to see that isolated galaxies sustain com- 
parable bar strengths. 

7. Summary and conclusions 

From the detailed morphological study of a sample of 45 iso- 
lated spiral galaxies, in H„ emission and optical red continuum, 
and comparison with numerical simulations, we come to the fol- 
lowing results: 

1 . The percentage of bars (60% in the optical) shows that iso- 
lated galaxies are not preferentially barred or unbarred galax- 
ies. The histogram of relative tangential force strengths is 
quite similar to other samples of galaxies selected without 
regard to any isolation criterion. 

2. Frequently, we observed a phase shift between gas and stellar 
components: the H„ emission is always leading with respect 
to the bar in the Gunn r image (shift angle ~ 10°). 

3. We interpreted the various global H„ morphologies observed 
in terms of the secular evolution experienced by galaxies in 
isolation. The main H„ classes can be related to the bar evo- 
lution phases. The observed frequency of particular patterns 



brings constraints on the time spent in the various evolution 
phases. 

4. We chose to group the galaxies presenting similar features, 
we defined three main groups: group E (19 galaxies) gath- 
ers the galaxies showing a strong central peak in the H„ 
emission, H„ emission along the spiral arms but not in the 
bar. Group F (9 galaxies) is constituted by galaxies with less 
gas, which do not present any central emission knot in H„. 
Galaxies in group G (8 galaxies) show Hq, emission in the 
bar. 

5. Other groups with fewer galaxies were also defined, most 
noticeable is the group EG (3 galaxies) featuring galaxies 
having characteristics in between those defined by the groups 
E and G. Hence, we could interpret these features as different 
stages of an evolutive sequence: G — > EG — > E — > F. 

6. Numerical simulations showed a predicted frequency of the 
G phase higher than the E phase, in contradiction with our 
observations. We attribute this discrepancy to a failure of the 
star formation recipe, since we used the usual local Schmidt 
law for the star formation rate. The frequently observed phe- 
nomenon of star formation avoiding the bar, in spite of large 
gas density there, suggests that the star formation law should 
depend also on other factors, in particular the relative veloc- 
ity of the gas in the bar. 
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